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ABSTRACT 

An intriguing fact about cosmic gas clouds is that they all appear to have neutral (atomic) hydrogen column 
densities smaller than 10 22 cm" 2 . Observations of damped Lya (DLA) absorption systems further indicate that the 
maximum A//// decreases with increasing metallicity. It is generally assumed that this trend is due to a dust-induced 
selection bias: DLA systems with high Nhi and high metallicity contain so much dust that the background QSO 
becomes too dim to be included in optically selected surveys. Here, it is argued that this explanation may not be 
viable. Instead, it is proposed that conversion to molecular hydrogen determines the maximum H I column density. 
Molecular hydrogen forms on the surface of dust grains and is destroyed by photodissociation. Therefore, the 
molecular fraction correlates with both the dust content and, because of self-shielding, the total hydrogen column 
density, and anticorrelates with the intensity of the incident UV-radiation. It is shown that the first relation can 
account for the observed anticorrelation between the maximum Nhi and metallicity. 

Subject headings: galaxies: formation — galaxies: ISM — intergalactic medium — ISM: clouds — ISM: 
molecules — quasars: absorption lines 



1. INTRODUCTION 

Hydrogen is the most abundant element in the universe and 
consequently the total hydrogen content is one of the most basic 
characteristics of gas clouds. Hydrogen is most easily observed 
' in atomic, neutral form through the Lya (1216 A) or the 21cm 
transition, either in absorption or emission. Because the line-of- 
sight extent of a gas cloud is difficult to measure and because 
the transverse size is generally unknown in case of absorption 
studies, it is not the mass, but the column density that is the 
main observable quantity characterizing the total gas content of 
a cloud. Perhaps the most basic observation about the distri- 
bution of neutral hydrogen column densities is that clouds with 
Nhi > 10 22 cm" 2 (> 80 M Q pc~ 2 ) appear to be extremely rare, 
if they exist at all. 

A well-known, possible explanation for the observed cut-off 
in the distribution of absorption line column densities is se- 
lection bias: if the dust-to-gas ratio is roughly constant, then 
higher column density clouds contain more dust, leading to 
a stronger extinction of background sources. Thus, the pres- 
ence of dust may cause sight lines through high column den- 
sity clouds to be missing from magnitude limited surveys (e.g., 
Ostriker & Heisler 1984; Wright 1990; Fall & Pei 1993). If 
dust-bias is important, then there should be an anticorrelation 
between the maximum column density and the dust-to-gas ratio 
(and hence metallicity Z) of gas clouds. Damped Lya (DLA) 
systems, i.e., absorbers with Nhi > 2 x 10 20 cm" 2 , do indeed 
show such an anticorrelation: Boisse et al. (1998; see also 
Prantzos & Boissier 2000) found that observed DLA systems 
satisfy [Zn/H] + log(Nni) < 20.5 and interpreted this as a dust- 
induced selection effect. 

The presence of dust in DLA systems has been established 
using two independent types of observations. First, it has been 
shown that QSOs with DLA systems in the foreground appear 
redder than those without DLA systems in the foreground (Fall, 
Pei, & McMahon 1989). Second, the relative abundances of re- 
fractory elements, such as chromium and iron, are significantly 

'in some extreme environments, such as broad absorption line QSOs, hydrogi 
measurements (e.g., Mathur, Elvis, & Singh 1995). However, these observations c. 



lower than those of elements such as zinc, which are thought 
to be only lightly depleted on dust grains (e.g., Pettini et al. 
1997a). 

Although dust seems to be present in at least a subset of 
DLA systems, the idea that dust-bias can explain the absence 
of high Nhi, high Z systems faces at least three potential prob- 
lems. First, 21cm emission line studies of nearby galaxies do 
not reveal higher H I column densities than are found in absorp- 
tion line studies. Radial profiles of the neutral hydrogen sur- 
face density of disk galaxies generally show a maximum value 
of max(A/#/) <; 10 21 cm" 2 (<; 8 M Q pc" 2 ) and always less than 
10 22 cm" 2 (e.g., Cayatte et al. 1994; Rhee & van Albada 1996). 
Second, preliminary results from the survey of radio-selected 
QSOs by Ellison et al. (2001), which is free from the dust-bias 
that affects optically selected samples, do not provide evidence 
for a previously unrecognized population of Nhi > 10 21 cm" 2 
absorbers. Third, the metallicities and dust-to-gas ratios of 
DLA systems are typically more than an order of magnitude 
below the Galactic values (e.g., Pettini et al. 1997a, 1997b) 
and Prochaska & Wolfe (2001) have argued that the implied ex- 
tinction corrections are far too small to explain the obscuration 
threshold proposed by Boisse et al. 

In this letter, a simple physical explanation for the cut-off in 
the distribution of H I column densities is proposed. Although 
the total hydrogen column density of a gas cloud will always 
increase as its density increases, the same is not true for the 
neutral hydrogen column density. As the (column) density in- 
creases, the fraction of hydrogen in molecular form increases 
and eventually the H I column density will stop increasing. Fur- 
thermore, the cross section for DLA absorption provided by 
self-gravitating clouds with high molecular fractions is low be- 
cause such clouds are compact and unstable to star formation. 
Because hydrogen molecules are formed on the surface of dust 
grains, the molecular fraction correlates with the dust-to-gas 
ratio (and thus metallicity). Hence, the conversion of H I to 
H2 naturally explains the observed anticorrelation between the 

columns in excess of 10 22 cnT 2 have been inferred from soft X-ray absorption 
not distinguish between atomic and molecular hydrogen. 



2 



A physical limit on the H I column density 



maximum H I column density and metallicity. In the remainder 
of this letter it will be shown that this idea works quantitatively. 

Earlier, related work on the conditions required for the con- 
version of H I to H 2 includes Hollenbach, Werner, & Salpeter 
(1971), Federman, Glassgold, & Kwan (1979) and Franco & 
Cox (1986). 

2. METHOD 

In this section I will describe and discuss the method used 
to derive the maximum neutral hydrogen column density as a 
function of metallicity for a given dust-to-metals ratio and inci- 
dent radiation field. 

Consider a sight line through a gas cloud of arbitrary shape. 
Let tin and L be the characteristic density and size of the ab- 
sorber. For now, let us assume that the cloud is self-gravitating, 
i.e., the pressure of the medium external to the cloud is low 
compared to its central pressure, and supported by thermal pres- 
sure. Such clouds will generally be close to local hydrostatic 
equilibrium, i.e., the characteristic size will be close to the lo- 
cal Jeans length Lj. If L ^> Lj, then the cloud will expand or 
evaporate and equilibrium will be restored on a sound crossing 
timescale. If L -C Lj, then the cloud is Jeans unstable and will 
fragment or shock to the virial temperature, equilibrium will 
then be restored on a dynamical timescale. Schaye (2001a) used 
this argument to derive the properties of Lya forest absorbers, 
to explain the shape of their column density distribution func- 
tion, and to compute their contribution to the cosmic baryon 
density. The results agreed very well with both observations 
and hydrodynamical simulations. Because the densities of in- 
terest here are much higher than those of Lya forest absorbers, 
the timescales for the restoration of hydrostatic equilibrium are 
much shorter and hence one would expect that the same argu- 
ment works at least as well for self-gravitating DLA systems as 
for Lya forest absorbers. 

In local hydrostatic equilibrium, the total hydrogen column 
density is given by the following expression (Schaye 2001a) 2 : 



(dP /drift > 0). Third, from these solutions the maximum Nhi is 
determined. The whole procedure is then repeated for different 
metallicities to derive the maximum neutral hydrogen column 
density as a function of the metallicity. 

Figure 1 illustrates the results for a model with metallicity 0.1 
solar and a dust-to-metals ratio of 0.5 times the Galactic value, 
values typical for DLA systems (e.g., Pettini et al. 1997ab; 
Vladilo 1998). The dust was assumed to have ISM properties 
and consists of a mixture of graphites and silicates (see the 
CLOUDY documentation for details). The gas clouds were ex- 
posed to the model for the UV/X-ray background radiation at 
z= 3 of Haardt & Madau (200 1) 5 , which includes contributions 
from QSOs and galaxies and yields H I and He II photoioniza- 
tion rates of 1.15 x 10" 12 s" 1 and 1.96 x 10" 14 s" 1 respectively. 
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where 7 = 5/3 is the ratio of specific heats for a monatomic gas, 
Y = 0.24 is the baryonic mass fraction in helium, f g is the frac- 
tion of mass in gas (excluding stars) and the other symbols have 
their usual meanings. To be conservative, f g will be set equal 
to unity. 

For a given metallicity, dust-to-metals ratio and incident ra- 
diation field, the neutral hydrogen column density is computed 
as a function of the density as follows. First, a grid of (NhAh) 
models is computed using the publicly available photoioniza- 
tion package CLOUDY 3 (version 94; see Ferland 2000 for de- 
tails), modeling the absorbers as slabs of constant density il- 
luminated from two sides 4 . The temperature is not a free pa- 
rameter. For each model the thermal equilibrium temperature, 
the mean molecular weight \i, and the neutral hydrogen col- 
umn density are computed self consistently. Second, the so- 
lutions (NH,nH,T(NH,riH),n(NH,n H )) are selected for which 
equation (1) is satisfied (i.e., Nh = Nh,j) and which are stable 

2 This expression differs by a factor 7T 1 / 2 from equation A5 of Schaye (2001a), who performed a purely dimensional analysis. Since this factor is usually included 
in the definition of the Jeans length, I include it here to obtain a more conservative upper limit on the maximum column density. 

3 http : / /www . pa . uky . edu/~gary / cloudy/ 

4 Illumination from two sides is approximated by doubling all column densities of a plane parallel cloud illuminated from one side. 

5 The data and a description of the input parameters can be found at http : / /pit to .mib . inf n . it /~haardt/ref model . html. The cloud is also 
exposed to the cosmic microwave background and to a cosmic ray density of 2 X 10~ 9 cm 4 , but these have no significant effect on the results. 
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FIG. 1 . — Hydrogen column densities as a function of density. The model as- 
sumes local hydrostatic and thermal equilibrium, metallicity Z = O.IZq, dust- 
to-metal ratio 0.5 times the Galactic value, and the model for the z = 3 UV/X- 
ray background radiation of Haardt & Madau (2001). There is a gap around 
n H ~ 10 -1 cirT 2 because there are no stable solutions in this region. 



Note that the curves in Fig. 1 do not show the evolution of an 
individual cloud. Each point along the curves corresponds to a 
cloud in local hydrostatic and thermal equilibrium, and higher 
densities correspond to lower cloud masses. If a cloud accretes 
matter, it will become Jeans unstable and fragment into pieces 
of higher Nh and lower mass which will again be close to the 
equilibrium curve. The model clouds become self-shielded at 
tin ~ 10~ 2 cm" 3 . There is a gap at nn ~ 10" 1 cm" 3 because 
clouds in this range are unstable (dP/dun < 0). At lower densi- 
ties the clouds are warm (T ~ 10 4 K) and extended (L ;> 10 3 pc), 
while clouds with higher densities are colder (r <; 10 3 K) and 

more compact (L <; 10 2 pc). For n H ^ 10 cm" 3 the clouds are 
highly molecular and clouds with higher densities have smaller 
H I column densities. The physical properties of self-shielded, 
self-gravitating gas clouds will be discussed in more detail else- 
where. 

The method described above assumes that the absorbers are 
self-gravitating, an assumption that may not hold for DLA sys- 
tems. Indeed, it is well known that the interstellar medium of 
the Galaxy has a multiphase structure and that H I clouds are 
generally pressure confined, although molecular clouds are not. 
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FIG. 2. — Metallicity as a function of neutral hydrogen column density. Model curves indicate the maximum Nhi (solid lines) and the HI column density for 
which the molecular fraction is 10 percent (dashed lines). Panel (a): From top to bottom the curves correspond to a dust-to-metals ratio of 0.1, 0.5 and 1.0 times 
the Galactic value respectively and the z = 3 UV/X-ray background of Haardt & Madau (2001). Panel (b): Curves are for a dust-to-metals ratio of half the Galactic 
value and the UV/X-ray background used in panel (a) multiplied by, from top to bottom, factors of 3, 1, and 1/3 respectively. Data points are the observed z > 1 
DLA systems taken from Pettini et al. (1997b; 2000), Prochaska & Wolfe (1999), de la Varga et al. (2000), and Molaro et al. (2000). The 1-ct errors are typically 
about 0.1 dex for both Nhi and [Zn/H]. 



Furthermore, it was recently argued by Schaye (2001a) that at 
least a subset of DLA systems appear to arise in galactic winds, 
in which case (ram) pressure is likely to be more important 
than self-gravity. Fortunately, the possibility that some (or even 
most) DLA systems are confined by external pressure does not 
affect the current analysis because a gas cloud with L < Lj will 
always have a smaller column density than a gas cloud in hydro- 
static equilibrium with the same density. Pressure-confined so- 
lutions will therefore lie below the solid curve in Fig. 1 and the 
assumption of local hydrostatic equilibrium is thus conservative 
when estimating the maximum Nhi- Recall that self-gravitating 
clouds containing stars and/or dark matter also have lower col- 
umn densities. 

Higher neutral column densities are possible for clouds that 
are self-gravitating and for which turbulent or magnetic pres- 
sure dominates over thermal pressure, and for clouds that are 
rotationally supported and oriented such that our line of sight 
is nearly perpendicular to their spin axes. One would, how- 
ever, expect these caveats to be important mostly for clouds in 
regions of ongoing star formation. In such regions the ISM is 
likely to be multiphase and the H I clouds pressure-confined, 
thus reducing their column densities compared to clouds in hy- 
drostatic equilibrium. Finally, it should be noted that although 
the integrated column density can be very large for a sight line 
through an edge-on disk galaxy, this column will generally be 
the sum of the column densities of a large number of clouds. 

3. RESULTS 

Figure 2a shows the observed [Zn/H] as a function of Nhi 
for DLA systems with redshift z > 1. The absence of points 
in the upper right part of the diagram, indicating a metallicity- 
dependent Nhi cut-off, was first noted by Boisse et al. (1998) 
who interpreted it as a dust-induced selection bias. The lack of 
data points in the bottom left part of the figure is most likely 
a selection effect: the corresponding zinc lines are too weak 
to be detected. Indeed, while higher quality observations of 
[Si/H] and [Fe/H] confirm the absence of high Nhi, high Z sys- 
tems, they do not show a deficiency of low Nhi, low Z systems 



(Prochaska & Wolfe 2001). 

The solid curves indicate the neutral hydrogen column den- 
sity in local hydrostatic and thermal equilibrium, computed us- 
ing the method described in section ^. From top to bottom the 
curves correspond to dust-to-metals ratios of 0.1, 0.5 and 1.0 
times the Galactic value respectively. Equilibrium solutions for 
pressure-confined gas clouds and self-gravitating clouds that 
contain a gravitationally significant amount of stars or dark 
matter are all located to the left of these curves. From the figure 
it is clear that the models can naturally account for the absence 
of high Nhi, high Z systems. A higher metallicity (dust-to-gas 
ratio) results in an increased Hi formation rate and, because of 
the increased cooling via metals and H2, a lower temperature. 
Both the higher H2 formation rate and, via equation 1 , the lower 
temperature contribute to the decrease in the maximum Nhi- 

It is important to note that the maximum column densi- 
ties are derived for single clouds, i.e., regions along the sight 
line over which the density is of the same order as the near- 
est (local) maximum, whereas the data points correspond to 
the integrated column densities of systems of clouds. Since 
essentially all DLA systems show substructure in their low- 
ionization metal lines, which are thought to trace the neutral 
hydrogen density, DLA systems must generally consist of col- 
lections of clouds. Unfortunately, the strength of the DLA line 
prohibits us from measuring the H I columns of the individual 
components. Hence, each observed Nhi value must be consid- 
ered as an upper limit to the column density of the dominant 
component. If, for example, the strongest component would 
account for 50 percent of the neutral hydrogen column (and if 
all components had roughly the same metallicity), then the data 
point should be shifted by -0.3 dex along the Afe-axis, which 
would place it comfortably to the left of most models. 

Although the models clearly show that there should be an 
anticorrelation between metallicity and the maximum Nhi, the 
exact Nhi values are somewhat uncertain. First, although self- 
gravitating clouds will generally not be far from local hydro- 
static equilibrium, there is no reason why they should be in ex- 
act equilibrium. Second, the derived column densities are sen- 
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sitive to the assumed intensity of the incident ionizing radiation. 
Figure 2b illustrates the effect of changing the amplitude of the 
UV/X-ray background. From top to bottom the curves corre- 
spond to models that use the same background radiation as was 
used in panel (A), but with the amplitudes multiplied by factors 
of 3, 1, and 1/3 respectively. The maximum Nhi increases if the 
background is stronger. Although these models are indicative 
of the uncertainty in the mean background radiation, the ioniz- 
ing radiation could be stronger around DLA systems if they are 
located near regions of star formation or QSOs. In such sys- 
tems higher Nhi values would be possible, although it should 
be noted that absorbers in these environments would likely be 
pressure confined. 

Extrapolation of the solid curves in Fig. 2 suggests that large 
Nhi values are possible at zero metallicity. The reason is that 
for the low-Z models the Nhi(iih) curves flatten off very slowly 
with increasing tin- The maximum Nhi is only reached when 
the density is higher than is typical of molecular clouds and 
when the characteristic mass of the cloud is lower than is typi- 
cal of stars. The models used here are probably inadequate for 
systems which are nearly fully molecular. In any case, such 
systems are so compact that the cross section for interception 
by a random sight line is very small and they are highly un- 
stable to star formation. Note also that the extreme low metal- 
licity models are not relevant for DLA systems which all have 
Z ^> 10~ 3 Z©. Nevertheless, it is instructive to look at a more 
robust quantity than the maximum Nhi, such as the Nhi corre- 
sponding to a fixed molecular fraction. The dashed curves in 
Fig. 2 show the Nhi at which the molecular fraction reaches 10 
percent, much higher than is typical of DLA systems (e.g., Pe- 
titjean, Srianand, & Ledoux 2001). These curves do turn over 
at low Z. For example, logNni = 21 .8 at zero metallicity for the 
ionizing background used in Fig. 2a. Because the molecular 
fraction rises rather sharply to /(H2) ~ 10" 1 (see Fig. 1) due to 
a thermal instability, the curves corresponding to /(H2) ~ 10~ 5 
(not plotted) are similar to the dashed curves in Fig. 2. 

Finally, the models cannot explain the apparent lack of Nhi < 
10 21 cm" 2 systems withZ > Z Q . If DLA systems with such high 
metallicities do occur in nature and if their cross sections and 
lifetimes are non-negligible, then it may be that a dust-induced 
selection bias still needs to be invoked to account for the ab- 
sence of such systems from current samples. 

4. CONCLUSIONS 

One of the most basic observational findings is that all gas 
clouds seem to have neutral hydrogen column densities smaller 



than 10 atoms per square centimeter. Studies of DLA sys- 
tems have revealed that the maximum Nhi decreases with in- 
creasing metallicity. It was argued that dust-induced selection 
bias may not be a viable explanation for these observations. In- 
stead, it was proposed that clouds with Nhi > 10 22 cm" 2 do not 
occur because the clouds turn molecular before reaching such 
high column densities. Furthermore, because clouds with high 
molecular fractions are much more compact and short-lived 
than clouds with low molecular fractions, the latter provide a 
much larger cross section per unit mass for DLA absorption. 
The maximum H I column density is a decreasing function of 
metallicity, mainly because the formation rate of molecular hy- 
drogen increases with the dust content of the system. 

It was shown that models of self-gravitating gas clouds in lo- 
cal hydrostatic and thermal equilibrium, with a dust-to-metals 
ratio somewhat lower than the Galactic value, exposed to a 
model of the z ~ 3 UV/X-ray background radiation, can roughly 
account for the observed anticorrelation between the maximum 
Nhi value and metallicity. The assumption of pure gas clouds in 
hydrostatic equilibrium is conservative because self-gravitating 
clouds containing stars and/or dark matter, as well as clouds 
that are confined by external pressure, have lower column den- 
sities. The model works particularly well if one takes into ac- 
count that observed DLA systems consist of multiple compo- 
nents, and that it is the integrated H I column density that is 
measured, whereas the model predicts the maximum Nhi f° r 
single gas clouds. 

Finally, the models presented here predict that the molecular 
fraction increases with increasing Nh and dust-to-gas ratio (and 
thus metallicity if the dust-to-metals ratio is roughly constant), 
but is anticorrelated with the intensity of the incident UV radia- 
tion. These correlations could be uncovered observationally if a 
sample is constructed in which two out of the three parameters 
(Nh,Z,Iuv) are roughly constant. There could, however, still 
be considerable scatter in these relations if only a fraction of 
clouds are purely gaseous and self gravitating. The models also 
predict a strong anticorrelation between molecular fraction and 
temperature that is fairly insensitive to the values of the other 
parameters. 
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